Universal aging properties at a disordered critical point 
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We investigate, analytically near the dimension due ~ 4 and numerically in d = 3, the non 
equilibrium relaxational dynamics of the randomly diluted Ising model at criticality. Using the 
Exact Renormalization Group Method to one loop, we compute the two times t, correlation 
function and Fluctuation Dissipation Ratio (FDR) for any Fourier mode of the order parameter, 
of finite wave vector q. In the large time separation limit, the FDR is found to reach a non 
trivial value X°° independently of (small) q and coincide with the FDR associated to the the total 
magnetization obtained previously. Explicit calculations in real space show that the FDR associated 
to the local magnetization converges, in the asymptotic limit, to this same value X°° . Through a 
Monte Carlo simulation, we compute the autocorrelation function in three dimensions, for different 
values of the dilution fraction p at Tc{p). Taking properly into account the corrections to scaling, 
we find, according to the Renormalization Group predictions, that the autocorrelation exponent 
Xc is independent on p. The analysis is complemented by a study of the non equilibrium critical 
dynamics following a quench from a completely ordered state. 
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The study of relaxational dynamics following a quench 
at a pure critical point has attracted much attention 
these last few years 0, S S 0- Although simpler to 
study than glasses, critical dynamics display interesting 
non equilibrium features such as aging, commonly ob- 
served in more complex disordered or glassy phases ^] . In 
this context, the computation of two times t, tw response 
and correlation functions with associated universal expo- 
nents has been the subject of numerous analytical as well 
as numerical studies 

In addition, it has been proposed ^ that a non trivial 
Fluctuation Dissipation Ratio (FDR) X, originally in- 
troduced in the Mean Field approach to glassy systems, 
which generalizes the Fluctuation Dissipation Theorem 
(FDT) to non equilibrium situations, is a new universal 
quantity associated to these critical points. As such, it 
has been computed using the powerful tools of RG, e.g. 
for pure 0{N) model at criticality in the vicinity of the 
upp_er _a:itical dimension due ~ 4 and for various dynam- 
ics 



An important question related to the physical inter- 
pretation of X in terms of an effective temperature |p 
Teff = T/X is its dependence on the observables [Tollllj|. 
In this respect, a heuristic argument suggests that, 
for a wide class of critical systems, the local FDR asso- 
ciated to correlation and corresponding response of the 
local magnetization should be identical, in the large time 
separation limit, to the FDR for the total magnetization, 
i.e. for the Fourier mode q = 0. This argument relies 
strongly on the hypothesis that the time decay of the 
response function of the Fourier mode q is characterized 
by a single time scale Tq ^ q^^, with z the dynamical 
exponent. 

Characterizing the effects of quenched disorder on crit- 
ical dynamics is a complicated task and indeed the ques- 
tion of how quenched randomness modifies these prop- 
erties has been much less studied. In particular, in this 



context of critical disordered systems, the question of uni- 
versality, i.e. the dependence of the critical exponents on 
the strength of the disorder, is a controversial issue [T^ . 
In this paper, we address these questions on the randomly 
diluted Ising model: 
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where Si are Ising spins on a d-dimensional hypercubic 
lattice and pi = 1 with probability p and with proba- 
bility 1 — p. For the experimentally relevant case of di- 
mension d = 3 fl^ . for which the specific heat exponent 
of the pure model is positive, the disorder is expected, 
according to Harris criterion to modify the univer- 
sality class of the transition. For 1 — p <C 1, the large 
scale properties of at criticality are then described by 
the following 0(1) model with a random mass term, the 
so-called Random Ising Model (RIM): 



where (p 



(p{x) and ip{x) is a gaussian random vari- 
able ip{x)ip{x') = Ad'^{x — x') and ro, the bare mass, 
is adjusted so that the renormalized one is zero. The 
static critical properties of this model have been inten- 
sively studied [13 both analytically, mainly using RG, 
within various schemes, and numerically [la |. The (per- 
turbative) RG calculations below the upper critical di- 
mension in d = 4 — e, which we will focus on here, con- 
firm the qualitative Harris criterion and predict that the 
critical properties of these models Q for different val- 
ues of p close to one, are described by a new disordered 
fixed point, which is independent of p. Therefore, an im- 
portant statement of this RG analysis is that the critical 
exponents, which can be computed in an expansion of 
■y/e, 6.17. 77 = 0(e), are universal, i.e. independent of p. 
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This was recently confirmed by Monte Carlo simulation 
in d = 3|l6j. over a wide range of concentration p above 
the percolation threshold Pc — 0.31. Although quantita- 
tive discrepancies were found with pcrturbative RG cal- 
culations [17i |. universality was demonstrated by taking 
carefully into account the (strong) corrections to scaling 
|16|. In the equilibrium dynamics, at variance with the 
pure case, the pcrturbative expansion of the dynamical 
exponent z differs from its high temperature value of 2 

already at one loop z — 2 = \J'^ + 0{e) indepen- 
dently o f jp, corrections to two loops have been computed 
in Ref. [13 and up to three loops in Ref. [201 . After a 
long debate, a recent numerical simulation '2 1| where cor- 
rections to (dynamical) scaling were taken into account, 
has confirmed the universality of z in d = 3, leading to 
z = 2.62(7) independently on the spin concentration p 
above pc- 

By contrast, much less is known about the non equi- 
librium dynamics of this disordered system at critical- 
ity. The critical initial slip exponent 0' vanishes to one 
loop I22I ■ and correction to two-loops have been com- 
puted |23|. This exponent has been recently computed 
up to two loops for the case of extended defects in [23 |. 
The two times response in Fourier space, "R-lt^ |23 and 
the correlation Cl^^ 12^, including the associated scal- 
ing functions, are known up to one loop. But although 
the dynamical RG predicts a universal value of the auto- 
correlation exponent Xc = d — z6' , for p close to 1, this 
statement remains an open question for a wider range of 
values of p. Furthermore, a non trivial FDR [2^, only 
for the total magnetization, was recently obtained to one 
loop, and it was argued, using the same aforementioned 
heuristic argument 0, to coincide with the local FDR. 
However, it was already noticed in Ref.j^^] that, due to 
the disorder, 'R■1^^ decays as a power law for q^t ^ 1. 
Therefore, the argument of Ref. is challenged for this 
disordered critical point, and, already at one loop order, 
the computation of the FDR needs a closer inspection, 
including an extension of the analysis of Ref. beyond 
the "diffusive" q — Q mode. 

In this paper, using RG to one loop, we obtain, for 
any finite Fourier mode g, the correlation function C^i^ 
and the FDR X^^^ , which are both characterized by scal- 
ing functions of the variables q^{t — t^) and t/tw In the 
asymptotic large time separation regime t^tw, the FDR 
reaches a non trivial value X°° , independently of (small) 
q. In addition, we explicitly compute the local FDR, 
which is a function of t/tw and reaches the same non- 
trivial limit Ar°° when t ^ t^^, which thus establishes on 
firmer ground the heuristic argument of Ref. 7] for the 
present disordered case. Besides, we perform a Monte 
Carlo simulation of the non-equilibrium relaxation of 
following a quench from high temperature with initial 
magnetization toq = at T'c(p) and compute the auto- 
correlation function. In the asymptotic regime, it takes 
a scaling form compatible with the RG calculations. By 
taking into account corrections to scaling, we show that 



the exponent Ac is independent of p. Finally, we compute 
numerically the autocorrelation function for the critical 
dynamics following a quench for a completely ordered ini- 
tial condition with mo = 1. We observe that the system 
is also aging and show that the decaying exponent Ac is 
strongly affected by this initial condition. 

We study the relaxational dynamics of the Randomly 
diluted Ising Model in dimension d = A — e described by 
a Langevin equation: 



ax,t) 



(3) 



where where {({x,t)) — and {C,{x,t)C,{x' ,t')) = 
2riT6{x — x')d{t — t') is the thermal noise and rj the fric- 
tion coefficient. At initial time ti = 0, the system is in 
a random initial configuration with zero magnetization 
mo = distributed according to a Gaussian with short 
range correlations 



[V?(x, t = 0)^{x', t - 0)], - T-H''{x - x') 



(4) 



Notice that it has been shown that Tq is here irrelevant 



(in the RG sense) in the large time regime studied here 
|26|. We will focus on the correlation Cl^^ in Fourier 
space and the autocorrelation Cu^ 
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{viq,t)ip{-q,t^)) , C'tt^ = {ip{x,t)ip{x,t^)){b) 



and the response TZft^ to a small external field f{—q,tw) 
as well as on the local response function R^^ respectively 
defined, for t > 



q ^ S{ip{q,t)) 
"™ Sf{-q,t^) 



^ S{ip{x,t)) 
Sf{x,t^) 



(6) 



where — and (..) denote respectively averages w.r.t. disor- 
der and thermal fluctuations. We focus also on the FDR 
ATj'^ associated to the observable ip l^: 
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L((^ — 1 at equilibrium. Notice also 
that for this choice of initial conditions Q), connected 
and non connected correlations do coincide for large sys- 
tem size. 

A convenient way to study the Langevin dynamics de- 
fined by Eq. Q is to use the Martin-Siggia-Rose gen- 
erating functional. Using the Ito prescription, it can 
be readily averaged over the disorder. The correlations 
(O and response © are then obtained from a dynami- 
cal (disorder averaged) generating functional or, equiva- 
Icntly, as functional derivatives of the corresponding dy- 
namical effective action F. This functional can be per- 
turbatively computed using the Exact RG equation 
associated to the multi-local operators expansion intro- 
duced in [2^ (2^ . It allows to handle arbitrary cut off 
functions c((7^/2Ag) and check universality, independence 



3 



w.r.t. c{x) and the ultraviolet scale Aq. It describes the 
evolution of F when an additional infrared cut off A; is 
lowered from Ao to its final value A; — > where a fixed 
point of order 0{y/e) is reached. In this limit, one obtains 
(for t> tw) from 



Kt^ and Ct^ 



dhl^tuHl,^^ = (8) 



dti 



(9) 

with fi(t) = — dtiTitti and where the self energy St^tj 
and the noise-disorder kernel i^titj are directly obtained 
from r at the fixed point. One finds: 



^"'^"Wi /j^"^*"*'^^' ^^^^ 

{la{t-t')-Jait + t')) (11) 



where 7a (x) — {x + a/(2AQ))^^. For concrete calcula- 
tions, we have used the decomposition of the cut-off fun- 
tion c{x) = J dac{a)e~'^^ = f^e'"^^. 

The computation of the correlation function C^^ re- 
quires the knowledge of the response, which we first focus 
on. By solving perturbatively to order 0{^/e) the differ- 
ential equation (|SJ|, similarly to what is done in Ref. [2^ . 
one recovers, in the limit g/Ag ^ 1 keeping the scaling 
variables v — q'^{t — t^) and u = t/tw fixed, the solution 
obtained in Ref. [23|, consistent with the scaling form 



-2+z+n 
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FR{q'{t~t^),t/t^) (12) 



where ^ = 5 



- 0{e) and the universal 30] scaling 

fimction Fr{v^u) admits also an expansion in powers of 

with 22] Fr{v,u) = F^^iv) = e-" \^(iv - 

1) Ei (v)e~'" + e~" — 1) where Ei(i;) is the exponential in- 
tegral function. At variance with the pure model at one 
loop 7], the large v behavior of F'^{v) is a power law, 
Fr{v) oc w"^, which already indicates that the heuristic 
argument of Ref.0 can not be applied here. Besides, 
when computing the local response R^^^ , one is left with 
an integral over momentum which is logarithmically di- 
vergent, indicating that this integral has to be handled 
with care to obtain the correct result, as the scaling form 
in H12(l is valid only for q/ Aq ^ 1. We thus solve pertur- 
batively the equation ||SJ| for any fixed q and obtain an 
expression for TZj^ consistent with the scaling form 



= ~9i{q) (^^j Fn{~g2iq)q^{t ~ t^).t/t^) (13) 

where ^1(9), 52 (9) are non- universal functions, i.e. which 
depend explicitly on the cut-off function c{x) and Aq, 



with the universal |30j small q behavior: 



dhTZl^Dt.t.nl 
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52(9) - 
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which thus allows to recover the previous expression in 
the asymptotic limit q/Ao <C 1 (|12(l . By computing the 
Fourier transform of Ti-u^ as given in Eq. I|13|l . we ex- 
plicitly check that the local response Rtt^, is consistent 
with the scaling form 
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^TCln(t 



) ( ^ 



2 (t- i^)l + ('i-2+'')/^ V*' 



(15) 



with Kd = Sd/i'^T^Y and where the non-universality is 
left in the amplitudes and A}j^: 



4" 



PR 



6e 
53 

In 



PR 



2A 



^R 



(16) 



At the order of our calculations ©(^/e), although z ^ 
2, this scaling form p5ll is compatible with local scale 
invariance arguments 31]. Notice also that, at this order, 
the scaling form obtained for R^^ could be written as 



(17) 



with a {d — 2 + ri)/z. Although this scaling form H17(l 
can not be ruled out at this stage, which would in princi- 
ple require a 2-loop calculations, it seems rather unlikely 
given the scaling form obtained in Fourier space l|12|l . 
where instead a logarithmic correction as in (|15|l is sug- 
gested by the large argument behavior of the function 

Fr»- 

We now turn to the computation of the correlation 
function in Fourier space Cff^ , which was only computed 
for the q = mode "25]. Solving one obtains an 
explicit expression, which, in the aforementioned scaling 
limit, is compatible with the scaling form 



a_ = T,q-^+'^ [ — ] Fc{q^t - t^),t/t^) (18) 



with the full expression 

Fc[v, u) = F^[v, u) + V^Fhiv, u) + 0{e) 
i;^0(i;,u) = e-^-e-^("+i)/("-i» 

Fh{v,u) = Fl?\v)-Fl?'^[v^^ 

2v \ , / 2v 



(19) 
(20) 
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where is the Euler constant. In the limit q ^0, our 
full expression (|19|l gives back the result of Ref. [23 . In 
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the large time separation limit m ^ 1, keeping v fixed, 
one obtains the result: 



Fc{v,u) = -Fc,M+0{u-^) 
u 



FcUv) = AcoovF^^iv) , Aco 



(22) 



(23) 



This ©(u"^) decay in H22II is expected from RG argu- 
ments, and has been explicitly checked for different pure 
critical systems Q. However, this relation H23() is a priori 
non trivial and can not be obtained from general argu- 
ments. This identity was also found for the pure 0{N) 
model at criticality to one loop order as well as in 
the glass phase of the Sine Gordon model with random 
phase shifts [s^l and it would be interesting to investigate 
whether such a behavior H23|) can be obtained from more 
general arguments. 

The full expression for C'j ^ (|19|l also allows to compute 
the structure factor Cli- It is obtained from H19|l in the 
limit u ^ 0, u — > 1 keeping v/{u — 1) = q'^t fixed, and 
we check that one recovers the previous result obtained 
in Ref. ;22]. Thus, one explicitly checks, at order ©(-^e), 
that the dynamical exponent z associated to dynamical 
equilibrium fluctuations is the same as the one associated 
to nonequilibrium relaxation. 

As noticed previously for the response function, the 
large v behavior of Fc{v,u) is a power law Fc(y,u) oc 
. Therefore, given the scaling form H18|l . the computa- 
tion of the autocorrelation Ctt^ has to be done carefully. 
Just as for the response, we thus compute the correlation 
function Cl^^^ from for any fixed q and then perform 
the Fourier transformation. One obtains the scaling form 



a, 



with 
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T{t/t^) (24) 
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The same remarks, concerning the response, made before 
Eq. (|17l) also hold here for the autocorrelation. 

We now turn to the FDR, first in Fourier space. Given 
the scaling forms for the response Ti-li^ H12|l and for the 
correlation Cli^ (|18|l we have explicitly checked here, the 
FDR X^^ takes the simple scaling form in the regime 
q/Ao < 1: 



the large time separation limit w ^ 1, keeping v fixed, 
one obtains, as a consequence of (P^ : 



0{e) 



(27) 



independently of v, i.e. of (small) wave vector q, which 
coincides of course with the asymptotic value for the q = 
mode obtained in Ref. j2^ . We can check easily, using 
the result of Ref. 7] , that this property, independence on 
V on the asymptotic limit, holds also for the pure model 
at one loop, and it was also found in the glass phase of 
the Sine Gordon model with random phase shifts ^32j . 

As we saw previously, the large v power law behavior 
of the scaling function F^{v) prevents us to use the ar- 
gument of Ref. [3 for the present case. Therefore on 
computes directly the FDR for the local correlation and 
associated response X^f^'^. It is also characterized by a 
scaling function of t/t^, which can be simply written as 



Tx{u) 



'(U+1)2 ^ 



6e / u — 1 



53 V w + 1 



(28) 
-0(e)(29) 



where Txiu) is a monotonic increasing function of u, it 
interpolates between 1, in the quasi-equilibrium regime 
for u — » 1, and its asymptotic value for m ^ oo given by 



lim (X- 



, _ , = lim (X,t°)-^=2+ j| + 0(6) 

(30) 

which shows explicitly, at order ©(Ve) that the asymp- 
totic FDR for both the total and the local magnetization 
are indeed in the same. 

Let us next present results from our Monte Carlo sim- 
ulations of the relaxational dynamics of the randomly 
diluted Ising model in dimension d — which were 
done on L X L X L cubic lattices with periodic bound- 
ary conditions. We first focus on the following situation 
where the system is initially prepared in a random initial 
configuration with zero magnetization mo = 0. At each 
time step, the sites are then sequentially updated : 
for each site ?, the move Si ^ —Si is accepted or rejected 
according to Metropolis rule. If one gradually decreases p 
the fraction of magnetic sites will be reached below which 
the system no longer exhibits a transition to ferromag- 
netic order at any finite temperature. This happens at 
the percolation threshold, for which Tc{pc) = ll^. 
For different values of p > Pc, we compute the spin-spin 
auto-correlation function defined as 



(31) 



Fx{q^t-t^),t/t^) 



(26) 



We have obtained the complete expression for the scaling 
function Fx{v, u), which at variance with the pure 0{N) 
model at criticality is a function of both q^t and q^t' . In 



In the following we will also be interested in the con- 
nected correlation function C(t,tw) defined as 

C{t,t^) ^l^J2{Si{t) Siit^n)) - {s^{t)){s^{t^)) (32) 
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FIG. 1: Log-log plot of auto-correlation function Ctt^ vs. t — 
t^. Inset: Scaling plot of Ctt„ as function oit/t^. Here, the 
system is initially prepared in a random initial configuration 
with zero magnetization. 



In order to obtain better statistics, Cu^ (or C{t,tu,)) is 
averaged over a suitably chosen time window At around t, 
with Af <C t. All our data are obtained for a lattice linear 
size L — 100, as an average over 500 independent initial 
conditions and disorder configurations. We also produced 
data (not shown here) for the spatial correlation function, 
for the same system size, to ensure that our results are 
not influenced by finite size effects. 

Fig. n shows the auto-correlation function Ctt^ as a 
function oi t — tw for different values of the waiting 
time = 2^1, 2^, 2^ 2"^ and 2* at p = 0.8. One ob- 
serves a clear dependence on tw, which indicates a non- 
equilibrium dynamical regime. We have also checked that 
for this choice of initial conditions, C{t,t^) and C{t,t^) 
do coincide. The scaling form obtained from the RG 
analysis H24|l suggests, discarding the logarithmic correc- 
tion, to plot 4^+''^/^C(t„ as a function of t/t^. Taking 
the values i] = 0.0374 from Rei.^ and z = 2.62 from 
Ref . [2ll| . we see in the inset of Fig. ^that, for p = 0.8, 
one obtains a good collapse of the curves for different 
tw Notice that such scaling forms are also obtained in 
more complicated disordered systems like 3-dimensional 
spin-glasses ^3^. 

However, for different values of p, the best collapse, un- 
der this form (|24|l . would be obtained for a p-dependent 
exponent (1 + t])/z. Thus one would conclude that this 
exponent is non universal Nevertheless, it is known 
[2ll| that such p-dependence occur due to corrections to 
scaling. Therefore, to include them, we extend the scal- 
ing form I24|) as 



1 



(^/ ) 



-Gp {t/t 



.(33) 

with b — uJd/z, where uj^ corresponds to the biggest irrel- 
evant eigenvalue of the RG in the dynamics, which is a 
priori different from the leading corrections in the statics 
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FIG. 2: t[^'^'^^^"Cu^/f{t - t^) as a function of t/U, for dif- 
ferent p = 0.5, 0.6, 0.65 and 0.8. Waiting times cor- 
responding to p = 0.5, 0.6, 0.65 are 2®,.., 2^" whereas for 
p = 0.8, tya = 2'',..,2®. f{x) is defined in the text. Inset 
: non- universal amplitude D[p) as a function of p. Here, the 
system is initially prepared in a random initial configuration 
with zero magnetization. 



[2l|. Unfortunately, we do not have any information on 
the function Gp{x). We will thus propose the simplest 
hypothesis, Gp{x) = Fp{x). In Fig. [21 we show a plot 

of - t^), with f{x) = 1 - Dip)x-'' : 

this results in a reasonably good data collapse of the 
curves for different t^j, for p =0.5, 0.6, 0.65 and 0.8. 
For each value of p, this data collapse is obtained via 
the fitting of 3 parameters : the exponents 5, z and the 
amplitude D{p). We found quite stable value of the ex- 
ponents z = 2.6 ± 0.1 and b = 0.23 ± 0.02, which are 
both independent of p. Our value of z, together with 
ujd = 0.61 ± 0.06 are consistent with the value obtained 
by Parisi et. at |23|. All the p-dependence is thus con- 
tained in the non universal amplitude D{p), as shown in 
the inset of Fig. |21 According to our data, the correc- 
tions to scaling in the quasi-equilibrium regime vanish 
for p = 0.8, i.e. D{p — 0.8) ~ 0, in agreement with a 



previous numerical c omp utation of the equilibrium auto- 
correlation function [ij. Notice that this value p = 0.8 
is also known in the statics, to minimize the correc- 
tions to scaling. 

As shown on the log-log plot in Fig. |2 and consis- 
tently with the RG prediction Fp{t/t^) (O de- 
cays as a power law for t ^ t^. However, this plot in 
Fig.[21would suggest that the decaying exponent depends, 
namely decreases, with p. We expect instead that this p- 
dependence is again due to corrections to scaling j21| . 
Consistently with the corrections we introduced in the 
quasi-equilibrium part of Cu^ in Eq. H33|l . we propose 
the form 

Fp{x) = A(p)a;(i+"--^'=)/" (l -I- B{p)x''') (34) 
where we (reasonably) assume that the dynamical cor- 



6 



10 



5 



O 



0.1 



0.01 



5 

„ 4 

a , 

< 1 


-1 





A(p) ^ 




B(P) 


G — — e ® T 




0.5 0.6 0.7 O.i 
P 



10 

t/t., 



100 



o 




0.03 



10000 



FIG. 3: Universality of Ctt„ for p = 0.6, 0.65, 0.8. The func- 
tion g{x) is defined in the text. Inset: Non-universal ampli- 
tudes A{p), B{p) as functions of concentration p. Here, the 
system is initially prepared in a random initial configuration 
with zero magnetization. 



FIG. 4: Log-log plot of the correlation C{t,tm) as a function 
of t — tw ior p — 0.65. Inset: Log-log plot of the connected 
correlation C{t ~ t^) as a function of t — tm. In the inset, we 
use the same symbols as in the main figure. Here mo — 1. 



rections to scaling are characterized by the same, p- 
independent, exponent b — 0.23 ± 0.02 as obtained pre- 
viously H33|l . Therefore, for each value of p one has three 
parameters to fit : the exponent Xc/ z and the amplitudes 
A{p), B{p). We obtain a quite stable fit for the different 
values of p, with the p- independent value of the decaying 
exponent Xc/ z: 

— = L05± 0.03 (35) 

z 

all the p-dependence being contained in the non-universal 
amplitudes A{p), B{p) (see the inset in Fig|21). As shown 
in Figl^l the curves for different values of p (and dif- 
ferent tw) in Fig. 121 collapse on a master curve when 
we plot tl'+''^^'Ctt„/(/(t - t^)g{t/t^)), with g{x) = 
A{p){l + B{p)x~^), as a function of t/tw This fact sup- 
ports universality of the long-time non-equilibrium re- 
laxation in this model. Our value for the exponent Xc/ z, 
together with z — 2.6 ±0.1 gives for the initial slip expo- 
nent 6' — 0.1 ± 0.035, which is in rather good a gree ment 
with the two-loops RG result O'^i^^^^ = 0.0868^. Al- 
ternatively, this exponent could be measured by study- 
ing the initial stage of the relaxational dynamics starting 
from a non-zero magnetization : this is left for future 
investigations [33 |. 

Here also, one obtains that the corrections to scaling 
in Eq. H34|l vanish for p = 0.8. We notice that this result 
is in apparent contradiction with the previous analysis of 
the non-equilibrium relaxation in this model performed 
in Ref. 21], where the focus was on the non-connected 
susceptibility, a one-time quantity, which instead claimed 
a "perfect Hamiltonian" for p ~ 0.63. However, the sta- 
tistical precision of our data does not allow us to make 
a strong statement about this point, which certainly de- 
serves further investigations. 



So far, we have focused on the relaxational dynamics 
occurring after a quench from a completely disordered 
initial condition, with zero initial magnetization mg = 0, 
to Tc{p). But it is also interesting to study how these 
agi ng properties depend on the initial conditions |3^ l3^ 
|37| . We have therefore performed numerical simulations 
where the system is initially prepared in a completely 
ordered state: 

S'i(t = 0) = +1 , V occupied site i (36) 

such that the initial magnetization is mo = 1. The sys- 
tem is then quenched at i = to Tc{p) and evolves ac- 
cording to the same aforementioned dynamical rules. We 
also compute the autocorrelation function C(i, tw) as de- 
fined in Eq. H31|) . The result of this computation for 
p — 0.65 is shown on Fig. 4, where we plot C(i,tu)) as a 
function of t — tw, for different tu, = 2^, 2^, 2^. Here also, 
one observes a clear dependence on the waiting time t^, 
which indicates that the system is aging. Notice however 
that, at variance with the previous situation (Fig. 1), the 
correlation for a given t — t^, decreases as increases. 
In addition, at variance with the previous case toq = Oi 
the behaviors of the connected C{t, t^) and the non- 
connected C(t, t^) correlations are qualitatively different 
: this is shown in the inset of Fig. 01 when one observes 
that C{t,t-uj) decays indeed much faster [s^. This prop- 
erty could be relevant for the computation of the FDR 
in this situation. The quantitative analysis of the cor- 
relation function C{t,tw) is shown on Fig. |31 Indeed, 
the curves for different t„ can be plotted on a master 
curve if one plots, for different tw, tw^^'^^^Cu^/ f{t — tw) 
as a function oi t/tw which suggests that also in that 
case the correlation function can be written under the 
scaling form as in Eq. with Fp{x) = Gp{x). How- 
ever, as illustrated on Fig. [S] the behavior of C{t,tyj) is 
strongly affected by the initial condition, the decay be- 
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FIG. 5: ti,~^^''^Ctt^/ f(t — tn,) as a function of t/t^ for p = 
0.65 and the two different initial conditions considered here, 
mo = and mo = 1. f{x) defined in the text. The straight 
hne is a guide line for the eyes. 

ing much faster when the system is initially in a random 
configuration with mo = 0. More precisely, as suggested 
on Fig. [SI our data for mg = 1 are compatible with the 
following scaling form 

C 1 / D{p) ^ft]"^ 

, t>i^ (37) 

where f3, v are the standard equilibrium critical expo- 
nents and where we have used the hyperscaling relation 
(ijv = (d — 2 + 77)72. Thus, although we can not show 
it analytically for the present problem, we believe that 
in that case of a fully ordered initial condition (toq = 1)j 
although the system displays aging, the exponent Ac is 
completely determined by the equilibrium exponents: 

Ae = - (38) 

V 

This relation (|38|l can be understood by considering 
C(t,0). Indeed, for this particular initial condition (|36|) . 
one has C(i,0) ~ M{t), where M{t) is the global mag- 
netization at time t. Therefore, at large time, from stan- 
dard scaling argument C(t, 0) ~ which thus gives 



the relation (|38() . Notice that this relation (|38ll is also 
found in the context of pure critical point (36i . i37i | . 

To sum up, we have performed a rather detailed anal- 
ysis of the relaxational dynamics up to one loop of the 
randomly diluted Ising model in dimension d = 4— e. The 
computation of the correlation function C^^ , including 
its associated scaling function, allows us to show that 
the Fluctuation Dissipation Ratio reaches, in the large 
time separation limit, a non trivial value X°°, indepen- 
dently of small wave vector q. Although, due to the broad 
relaxation time spectrum induced by the disorder, the 
standard argument of Ref. can not be applied here, 
we have performed an explicit computation in real space 
which shows explicitly that the limiting FDR associated 
to the total magnetization, on the one hand, and the 
local one, on the other hand, do coincide. And in this 
respect, it would be interesting to further investigate the 
FDR associated to other observables, like the energy for 
instance . These properties could also be tested 

in numerical simulations. 

In addition, we have computed numerically, in c? = 3 
the autocorrelation function. It is characterized by a scal- 
ing form fully compatible with our one loop RG calcula- 
tion in real space. We have however shown that this two 
times quantity is strongly affected by corrections to scal- 
ing, which remain to be understood more deeply from an 
analytical point of view. By taking them properly into 
account, our data suggest a universal, i.e. p-independent 
autocorrelation exponent Ac, which provides an "indi- 
rect" measurement of the initial slip exponent 6' , which 
is in reasonably good agreement with two-loops RG pre- 
diction. Finally, we have shown that the critical dynam- 
ics following a quench from a completely ordered state 
(toq = 1) displays also aging, but with a quantitative 
different behavior, the decaying exponent Ac being in 
that case completely determined by the equilibrium ex- 
ponents. 
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